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Abstract 

Nonperturbative spin-dependent forces of quarks in a baryon are 
calculated directly from the QCD Lagrangian in the framework of 
the Field Correlator Method both for heavy and light quarks. Re- 
sulting forces contain terms of 5 different structures, only one being 
known before in asymptotic form.Perturbative terms obtained by the 
same method are standard and have different signs and structures with 
respect to the corresponding nonperturbative ones, implying possible 
cancellations for some baryonic states. 

1 Introduction 

The spin structure of baryons presents a still unsolved problem,both on par- 
tonic and quark model level. For excited baryon spectrum the apparent 
small spin-orbit splitting of some baryonic states is a topic of vivid discus- 
sions Some baryonic states, like Roper resonance A/" (1440) or A(1405) 
are not yet explained in the traditional framework of relativistic quark model 
(RQM) p|, IJ.More detailed information about the spin structure of baryons 
comes from the polarization experiments on electroproduction of excited res- 
onances, || || which effectively measure the convolution of the baryon wave 
function, and is very sensitive to its structure. 

Meanwhile the theoretical knowledge of the quark spin forces in the 
baryon is limited to the perturbative expressions calculated decades ago fTj, 
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and nonperturbative spin-orbit Thomas term, written in the framework of 
RQM f§. 

In applying these results to light baryons the notion of constituent quark 
masses is introduced in RQM,these appear in spin-dependent (SD) forces and 
play role of fitting parameters. 

It is the purpose of the present paper to derive SD forces in a baryon 
in a most straightforward way from the QCD Lagrangian with the nonper- 



turbative vacuum described by vacuum field correlators || [T(J- Limiting to 
the lowest (Gaussian) field correlators one can express all terms of SD forces 
through the scalar functions D and D\ representing this Gaussian correlator 



II]. High accuracy of such procedure is supported by recent lattice data [12 



and contribution of higher correlators can be estimated to be of the order of 
few percent |fT3| . 



The functions D, D\ are themselves measured on the lattice [H| and also 



found in analytic approaches |15|, |T( 

An essential element of the present approach is that it is not connected 
to the heavy mass expansion, and can be applied also to light quarks in a 
baryon. In this case an effective Hamiltonian is constructed from first princi- 
ples, which contains einbein (auxiliary) fields. It was shown previously, that 
stationary point of these einbein fields yields exactly the constituent quark 
masses which can be expressed unambigiously through the only parameter 
of this approach - the string tension a. 

The decisive check of this procedure is the calculation of baryon mag- 
netic moments, since they are inversely proportional to the quark constituent 
masses. That was done in JL7[] and results agreed with all known experimental 
data within ~ 10%. 

The SD forces derived below in the paper are computed as a series in field 
correlators (cumulants) with growing powers of fields. 

The lowest (Gaussian) term yields SD forces inversely proportional to the 
square of constituent quark masses. Having in mind the high accuracy of the 
Gaussian approximation [13| and baryon magnetic moments [17] one should 
expect that SD forces found below have the accuracy of the order of 10%. 



Analogous expressions for heavy quarkonia [18, 19 1 and light mesons [20 
have been reported earlier and estimated for realistic meson system respec- 
tively in m and H . 



The plan of the paper is as follows. In chapter 2 a general expression 
for the 3q baryon Green's function is introduced, and the Fock-Feynman- 
Schwinger (world-line) formalism is used to reveal the dependence on gauge 
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fields with spin operators explicitly written. Averaging over those with the 
help of the Field Correlator Method (FCM) one finally obtains an expression 
for the Green's function written in terms of field correlators. 

In chapter 3 a special case of heavy quark masses is considered and all SD 
forces are obtained in closed form, expressed in terms of correlator functions 
D and D x . 

In chapter 4 the perturbative contribution to SD forces is written down. 
In chapter 5 a general case is considered when current quark masses can also 
be vanishingly small, and SD forces are written again in terms of integrals 
over functions D and D\ with constituent (dynamical) masses entering in 
the denominator. 

Chapter 6 is devoted to the discussion of the relativistic structure of SD 
forces in the excited baryon spectrum. Comparison to other results in the 
literature is also made and possible extension of the method is suggested. 
Main points of the paper are summarized in Conclusions. 



2 3q Green's function with spin insertions 

Following |], ID, ^3| we consider the 3q Green's function, which can be written 



as 



G3g(x,y) = tr L 



Font II( m * 
i=l 



D) J Di2 i Dz^e- Kl (W 3 expgaF)T i 



(1) 



where tr^ is the trace over Dirac matrix indices, T out (r in ) are final (initial) 
state operators creating given 

jPC 

assignment to the 3q state and we have 

also denoted as in ]17|, ^ 

,2 



dt 



'"" +4V(t) + l) 



{W^exp(gaF)) = try exp 



E 

Ln=0 



nl 



((F(l)...F(n)))dp(l)..dp(n) 



(2) 



(3) 



Here and in what follows F(l) is always implied to be gauge-transported 
to one point x ,namely F(l) = F(z(l),x ) = $(x Q ,z(i))F(z(i))$(z(i),x ) 
where <&(x,y) is defined in the Appendix, and finally, dp{n) = X)f=i dp^(n), 
with 

1 u\ dt„ , . 

(4) 



2pi{t n )' 
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The integration in (Q) extends over all 3 lobes of the minimal area surface 
(Si + S 2 + S3) inside the quark trajectories z®(t) and the string-junction 
trajectory z^ Y \t). We have also denoted a^ u = j^il^lv — 7^7^) and every- 
where Euclidean space-time is used (till the last moment when the resulting 
Hamiltonian is obtained in Minkowski space-time) with 7-matrices 

74 = 7o = (3; 7i = -iP<*h lp.lv + Ivln = 26^. 
Note also that notation try means 

tr Y P = -eabc^a'b'c'Pabc/a'b'c'- (5) 

In the combination F(k)dp(k) in (|3|) one can write 

^2 = (^)e; "«b)> * =w - ^ 

As it was shown in fTO] the spin-independent part of (0) which obtains ne- 
glecting the S term in (f|), yields at large quark separations, — z^| ^> 
T g , i — 1, 2, 3, the familiar area-law asymptotics 

(W 3 ) = exp[-t7(5i + S 2 + 5 3 )] (7) 

implying linear confinement for each quark. In what follows we shall use the 
general expression (1) to derive the spin- dependent part of interaction both 
for heavy quarks (expansion in inverse powers of mass) and for light quarks. 



3 Spin-dependent interaction in 1/ra expan- 
sion 

To illustrate the method we shall start with the derivation of spin-dependent 
(SD) forces via 1/m expansion. Defining the SD potential as Vsd, one can 
write Gs q ~ e - TV so ^ \ _ TVsd, and for Vsd the following general form will 
be obtained below, similar (but not identical) to the corresponding form for 
heavy quarkonia pi, IS , 



B(2) w( 3 h ^-«L (i) / 1 d Vl , 1 de 



Vs D (r« r (2 U (3) ) = £ 



2m, 2 {RWdR® 2.R«dR« 
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+ N C - 1 f£ 2m imj W) dRW 

+ E 

Kj 

where n = ^j-,'R i j = Rj — Hj. We assume that current quark masses are 

large, y/a, i = 1,2,3, and hence also are large, since the latter 

are defined through rrii and a in the stationary point analysis JK| [L7| and 
always satisfy /ij = rrii + 0(l/rrii). Hence for simplicity we keep in the 
following fii = rrii ^> \f& and expand in inverse powers of 1/rrii. 

As it was observed in [18|, [19] the SD terms of the lowest order (1/mf , m m . ) 
come from 3 different sources: 

A) Diagonal terms in (||) are kept together with diagonal terms in A< = 
(mj + yielding one power of 1/mj. An additional power of ^- or ^- 

then comes from the expansion of (W3). This yields spin-orbit terms V[, 
and V5. 

B) The off-diagonal terms are kept both in (^) and in Aj. This gives 
spin-orbit potential 

C) Diagonal terms from two matrices (^J) with % 7^ j are retained. This 
yields spin-spin potentials V 3 and V4. We now calculate the SD contributions 
from A) - C) point by point. 

A) From fl3|), (f|) one gets for i = 1 

(trv[(l + ^ B k {z^\l^)dt 1 )W z \) « 1 - TF« (9) 

Mi jo 

Using the relation igF UI/ W = x —, , W which obtains easily with nonabelian 
Stokes representation for W, one has 

(tr Y F Xa (x, z )W(C)) = tr Y {ig J ds^F^z, z )F Xa (x, z )W 3 (C)} (10) 

one can rewrite the l.h.s. of @ as 

trrW+f-oFf dh(tr Y B k {z^\t^) [ ds^(u)F^(u,x )W 3 (C')). 

ZiTYIj\ J J S(Cj 

(11) 

In ([TO]) and (|TT|) the common reference point z is chosen to make both 
expressions gauge invariant; as it will be seen, this point will not appear in 
the final equations. 

5 



)V±{Ri 



3(ffWn)(ffWn 



12mimj(N c 



+ 



l2mimj{N c - 1) 



+ V 5 , (8) 



In ( |TTD the contour C is deformed due to orbital momentum of quarks 
as compared to the zeroth-order contour Cq consisting of straight lines. This 
is essential since otherwise the vacuum average (B k Wz{Co)) vanishes since 
it is odd with respect to reflection Z{ — » — Z{ i ^ k. Therefore all nonzero 
contribution in (11) is due to deflection of the quark path in C from the 
straight line in Cq. 

At this point we shall describe the quark trajectory z£'(t) and the cor- 



responding string piece Wjfi from the quark position to the string junction 
(which we for simplicity take at the origin). 

wf{t,P) = zf{t)P, 1 > /? > (12) 

= d^dte ik J^- (13) 
where L$ is the (Minkowskian) angular momentum of the z-th quark 

Lf=im i e am tffz®, R«=z«-z( y ) = z«. (14) 

Similarly da k l \ = R^dft^du^ and one arrives at the result 

c /9( i ')r( 1 ') 
(B k (z®,ti)W 3 (C)) = ig / dp^du,^—^(B k B n (u 4 ,(3)W) + 

J irrii 

+ igj d(3 {i) du A R ( f\u i ){B k E l {u i ,(3 {l) )W). (15) 

Denoting 

dA 

(S fc (z (1) ,ti)^(u,«4)W3) = e kin {u n - 4 1} )^ (16) 

one obtains 
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(tWl^ r T dAn 

-wl *»<*-Mh£ (17) 

where we have used the relation u n (u A ) — z$(t\) = z^{u a ){ua — ti). 

Until now we have not used the Gaussian dominance of the vacuum, i.e. 
the fact that (W 3 ) is saturated by the lowest cumulant (FF), which is found 
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to be an accurate approximation [IB]. Using it one can write (FFW 3 ) 



(FF)(W 3 ), and introduce scalar functions D,D X for tensor (FF), as it was 
done in [ |TT| . 

Referring the reader to the Appendix for the corresponding relations, one 
finally obtains 

fldVA® f R d\^ A, r°° , rr ^ /x . ^ , , 2 9D U 



RdRJ Jo R' R J— 00 v ' ' d\ 21 

J* dV Jo RM' (18) 
Till now we have taken into account the interaction of the spin of the i-th 
quark with the surface Si, which yields the term (^V/)w multiplied with 
(m?) _1 . At this point we consider the interaction of the z-th quark spin with 
the (deformed) surface Sj, which will give the term V{ in (^j). For this one 
needs to consider a vacuum average of two F's from two different surfaces Si 
and Sj. 

In general we have for two F's transported to the same point x (a, ...rj 
are fundamental color indices) 

(tr FF) 1 
(F(u, x) a ^F{y, x) rn ) = — _ - (S av S^ - —S^S 7V ). (19) 

Taking into account (|5|) and the relation 

tr Y <f>aa>(xiy)$pf3>(x,y)®w'(x,y) = 1 (20) 

one obtains 

,t- / i . (trF(u,x)F(v,x)) , , 

tr Y (F aa ,(u,x)Fw(v,x)) = K - ^ e(JVe _\j (21) 

where we have also accounted for different orientation of plaquettes in Si and 

Sj. 

Now proceeding as in flT5| ) one has 

(B^W 3 (C')) = g [ l P^dpWdu?^(B®B%Xu l £\p (J )R(j))w 3 )+ 
J Jo rrij 

+ ig J dl3 {j) du < i ) R\ 1 \u i i ) ){B l i ) E ( i l) W 3 ). (22) 



At this point one can use Gaussian dominance and relations ([H]) to obtain 
finally 

1 T a//t->\A f 1 jn(i)n(i) f°° j_.rw._fMl A , n , //(ii)\2 , 2\ ^1 1 



(^'(*)) . = I JZ dv[D(rM, u)+Dx + ((r«>) 2 + u 2 ) 

(23) 

In a similar way one obtains from the first term on the r.h.s. of (|22"D ,with 
the use of the last term on the r.h.s. of (Al). 

r (^Wl^lff, (24) 



Here we have defined = R^ — /3^R^, and one should take into 
account, that .D, L>i depend on their arguments as 



D(r,u) = D(Vr 2 + u 2 ). 

This concludes derivation of terms with the procedure A) and one goes 
over to the next point. 

B) Following (||) one can write for the corresponding term of a given 
quark (i) 

VsdT ~J2^ { { <r (i) P W rrn-in JV^E ) W >Wh ( 2 5) 
Now one can use relation 

(E k (z^^)Ws(C )) = . * (W f% ]) m = - 1 i(^(C )) (26) 

igda k4 {z^),t( 1 )) igdB> k > 

where the following notation was introduced for the spin-independent poten- 
tial e(R( x \R( 2 \ i? (3) ). 

(Ws(C )> = exp[-s(R^\R^,R^)T}. (27) 

Note that in (p5| - p7|) one can keep in W$(C) the unperturbed (straight- 



line) contours for quark trajectories since the pref actor in (25) is already 
0(l/m 2 ). Keeping in mind relation 



fEf^vf = E®p® + ie kln a®E®pP (28) 
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one recovers the second term on the r.h.s. of (^). 

C) Here one considers spin-spin interaction and the corresponding term 
looks like 



1- V S T>T = tr Y y E((l + ^- J <T«BW(^U)^)x 

i>j * 

x (1 + ^- J tr^B^iz^^dt^Ws) (29) 

Identifying spin-spin terms in and using and relations for (BBW3) 
in Appendix one arrives at 

•tAao) V"^ / aVCJ k a k> i x mini 1 \ 



4m i m i (iV c -l) L v v ' <9(u) 2/ 

where notations are used 

u = R w -R (i) , u = U-tj- (31) 

Rewriting ( |30|) as 

VsB ~h 12^(^-1) (32) 

V A {u)= dv(3D(u,v) + 30^11,1;) + 2u 2 —^) (33) 

J —00 <7M 



one has 



V 3 (u) = - I duu 2 £J \ (34) 



Sij = 3(<7 (i) n)(<x^n) - o^V^, n = A- (35) 

This concludes definition of all NP spin-dependent terms in (|8|) to the 
order 0(l/m 2 ) and in the approximation when only lowest, (FF), correlator 
is retained in the Wilson loop. 

Now comparing our expressions for V{, V^', e', V3, V4 with the correspond- 
ing ones for heavy QQ case, given in [IS, , one can see that they coincide 
exactly, the only difference being that one should sum up over all 3 quarks 
for V and e', and take a double sum, i < j, for V 2 ', V3, V4. In addition there 
is a term V5 which is of 3 body character and vanishes in two-body situation 
since in that situation L^'r^ = L^r = 0. 
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4 Perturbative spin-dependent forces 



We are now in position to consider also perturbative contributions to the SD 
potentials, which were calculated in @ . The easiest way for us is to remem- 
ber, that to the lowest order, 0(a s ), all perturbative terms are pair-wise 
interaction of quarks, and they can be reconstructed from the expressions 
obtained above, Eqs. ([18]), (|23|), (H),©-©, usin S °( a s) contribution to 
Di(x), while D(x) does not have contributions at this order, f25fl 



D(x) = D NP (x), D x {x) 



3^x~ 4 



NPi 



(36) 



It is rewarding to realize that D\ does not enter into V{ (terms containing 
Di in (|TBp cancel exactly), so that perturbative contribution occurs only in 
the nondiagonal, i ^ j, terms in (|8|) and in e' . 

In the QQ case analogous calculations have been done and compared to 



standard ones in |19 

We start with the ^ term and rewrite it in the original form (|26|), not 
assuming that e depends on only, but also on — , as it is for the 
Coulomb term to be added to e. In this way one obtains from 



SD 



4m? dR 



(0 



(37) 



and for e —>■ e + V c 



couh 



V, 



coul 



2a 3 



>j \R(V-RU)\ 



one obtains 



(e,peri) 



SD 



2a, 
3 



E 



X P 



www 



4m?(RWy 



+ 



x P 



0')W(j) 



4m?(i2«?); 



(38) 



This expression coincides with the corresponding in 0. Consider now 
nondiagonal spin-orbit term, equivalent to V^. Instead of using replacement 
( 36|) and doing integrations in (|23|), we start from more general expression 
(ISI) to derive 



v (2,pert) _ 
^SD 



2a« 



3(iV c 



E 



[<r® (R fe) x p 0')) + <r^(R w x pW)] 



rriim^RWy 



which is again in agreement with 



(39) 
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Next we consider the spin-spin interaction. Here it is straightforward to 
replace in V 4 © D 1 as in © D 1 -> Di pert) = gsj to obtain 

vt H \r) = f ^(3L>? er * } (r, + 2r 2 ^H = 
J-oo or z 

= ^%(3) (r) . (4Q) 




In a similar way one obtains for V^'^. 



Vt rt \r) = % (41) 



One can also persuade oneself that V5 has no pertubative counterpart 
since there — > and it is orthogonal to and therefore the total 
perturbative SD contribution to the order 0(a s ) can be written as 

.Avert) _ TA (e,pert) T/ (2,pert) \ - (T {t) (T^Vjf^ (R^) + SjjV^^ (R^) 
V SB - V SB + V SB +2^ io mm (AT T\ \ > 



12mi7rij(N c 



where explicit form of 4 terms on the r.h.s. of (|42l ) is given in (|38D , (|39l), ( |10| ) 
and (|I|). 

Our results for V^f ert \ V^ pert \ v}^^ coincide with the corresponding ex- 
pressions in m , however our V^ pert ^ is 2 times smaller than the corresponding 
term in |§. In the next sections we shall argue that for light quark this term 
gets indeed twice as big, since there one should replace m; — > and for 
the coefficient appears to be 2 times larger. 



5 Spin-dependent forces for light quarks 



In sections 3, 4 the SD forces have been obtained as an expansion in 1/rrii, 1/rrij 
taking all 3 quark current masses large, rrii ^> y/a, i = 1, 2, 3. 

It was noticed before [20| however, that general expressions ([]])- dl) f° r 
Green's functions written in FSR, with the einbein function Hi(t) introduced 
as in [10], allow to obtain expressions for SD forces also for light quarks 



without 1/m expansion, and the corresponding terms for the meson case 
have been written before [|I0|, |22|. Below we demonstrate in this section 
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that the same procedure works also for the 3q case with light current quarks 
masses as well. 

We start again with general form (|3]) and instead of expansion in 1/m 
(or which is equivalent for heavy quark) shall do the only approxima- 
tion: keeping in the sum in the exponent (|3]) the lowest (Gaussian) cumulant 
((F(1)F(2))). This approximation was recently supported by lattice data for 
Casimir scaling while higher cumulants provide (for Wilson loop) less 



than 2% p|. 



We start with the spin-spin interaction, which is easily obtained keeping 
in (|3p the bilocal term and in ([|) only the a- dependent term. In ([3]) one 
should take into account that F(l) and F(2) belong to different lopes Si, S 2 
of the S123 surface, hence n = 1 for each of them; moreover one uses (|H5p and 
(|2l|), the latter with opposite sign, since orientation of S(l,3) and S(2,3) is 
the same in our case. As a result one obtains 



The combination o-WcrW in © is a product of two 4x4 matrices, which can 
be split into the product of Pauli spin matrices <Tj and chiral 2x2 matrices 

1 and pi = [ ^ o ' ' ^ 1US one can rewr ^ e dH) as 



= £ r An (N \ °$°$[(*B m {f)B n mi >< 1) + 
i>3 z ¥nA f iV iV c — l ) 

+ (trE m {i)E n {j))(p 1 x Pl ) + (trB m (i)E n {]))(l x 
(trE m ^B n U))( Pl x i)] = V&\BB) + V^\EE) + 

V^\BE) + V^\EB). (44) 
Using formulas from Appendix for correlators of B, E one has 

where u = R w - R (j) . 



12 



One can see that coincides with (]3l3-|3l>|) with substitution 



m i , rrij — > 
Hi,IJLj- For EE term one obtains 

^SD ~ 1^ 19 7T7 77 (Pi X Pi) (4b) 

where we have defined 

V A (u) = / t&/(3D(ti, v) + 3D 1 (u, v) + (3z/ 2 + u 2 ) l K ' ; ) (47) 

Finally for the last two terms in fl44] ) one has 



( ct « xo -(;)U(P^ + p^) r oo 9jDlfu ^ 

^SD l-D-C'J - -^SD (-C'-dJ - "j 777 77 / TT2 v dv - 

(49) 

This concludes calculation of spin-spin interaction. 

We turn now to the calculation of spin-orbit terms. The corresponding 
expression in (0) can be written as 



{3 r dt- 
J2 i J ^il ds p ] A U ) D ^,pa(z,u) + 



w^i^.S ^ dsb M°>^M\ (50) 



In ( |50|) we have defined as in Appendix 

g 2 

D^ VtPa (z,u) = —(trFp U (z(ti))F pa (u)) = D(h)(5 w 5 ua - 5^5 up ) 
+~[dphp6 1/(T + perm.}Di(h) 

and hp = z^ti) - Up. 

The Dirac structure of the exponent in ( |50D is a sum, which can be written 
with notations from ( 44]) as 

y(so) = J2 V^diag)^, x J ^ + v {so,nondia 9 )^{{) ^ J . ^ ^ ^ 

i=l 
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In the 4x4 matrix cr^ we first consider the diagonal part. Repeating all 
the steps leading to (|18D and (p3|), ( [24]) one has the same expressions with 
the replacement — - > /ij , i = 1, 2, 3, namely 

v {sodiag) = j^ V^lP 1 1 <t(0lW + ^'l'" 1 dV 2 \ 

SD ]fe 2/i| RP>dRP> + N e -ifz 2 W mam] 

(52) 

Let us now consider the nondiagonal part in ([51]), which can be written 



as 



r 3 p ^.j- 

(W 3 exp(gaF)) {s0inondiag) = exp y ^^(D^^^dsfjlt/ ) 



(53) 

Now taking into account (0), (0) one has writing (|53"D in the form 



(W 3 exp(gaF))( SOtnondiaa ) = exp 



i 3 

-T </ V 4- V V {ij) 

1 / so,nondiag) ' / / (so,nondiag) 

i=l i<j=l 



(54) 

one finds from (f)3|), replacing D^ vpa from (Al) and using flHf), (]L4]) 



T/(") _ a(*») I I A(") 

^(.so.nondiag) — EE EE "+" Q EB l 00 J 

and similarly for V^ nondiag ) . For terms on the r.h.s. of (p5|) one obtains 

A (!i) = / dvdu(D(v, u) + D a + z/ 2 ^-) 
jo az/^ 

(<ru)(R ( ^^M^3# (i) , u = R (i) /?; 
an" 1 

= -\f ^ ^iKfr) A 1 ^^(lW x ^WjpW. (56) 
2/if J~oo Jo 
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Similarly for A^' one obtains from the second term on the r.h.s. of (53) 



. m . <rR (j) AW) 



2^i#) (JV C - 1) : 



/*1 /*oo 7~) 

JO J-oo OV l 



2/!j/!j J-oo CW Z JO /ij \lj 

(57) 

Here we have defined r^) = R w - ^R^\ 



6 Discussion 



Let us now discuss the results obtained in the paper. For the heavy- 
quark case the nonperturbative dependent potential is given in (§), and the 
perturbative part in (f22|), so that the total SD potential is 

-rrftotal) rAnonpert) T /(>ert) / ro\ 

V SD - V SD + ^SD • \p°) 

Perturbative part agrees with that obtained long ago in and repeated 
in many subsequent papers. However in |J the term V^ £,peTt ^ is taken twice 
as big as fl38f) (or corresponding term in 0). A possible modification for 
light quarks, which can produce this increase is discussed later. 

The nonperturbative part y r g ^ on P ert ) (g) consists of six terms, which were 
never fully written before. 

Only asymptotics at large distances of the first term in (ffl) has been 



written before in 271 and later in [261, one can find it from 



V (non P ert ){Rii) ^ oo) = _^^_ (59) 

In m were instead postulated the pairwise nonperturbative spin-orbit forces, 



which contradict expressions derived in this paper and in |26[. All other 
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terms, proportional to V^', V 3 , V4 and V 5 have never been written for the 3q 
case, while for the qq case the corresponding terms (except for V5) have been 
written in ]I8 |, ||19|| . The term V5 which has no counterpart in the qq case, is 



completely new, and its physical implication is still unclear. 

We now turn to be light quark case. Here the total SD "potential" is in 
general a sum of product of (4 x 4) matrices, which can be written as 

{Might quarks) _ V^/T/(ij)i v i i f/fe) \ 

V SD — Z^\ V diag L i X L j "+" V nondiag.) l DU J 

where Vnondiag contains terms like 1 j x , p u x 1 ^ , x p X j > an d Pi 

where each entry in pi, is 2 X 2 unit matrix. 
Now for vjfag one has 

V^ = V&\BB)+V&**> (61) 
where the first term on r.h.s. of (|H]) is given in (^) and the second in 

m 




(|52|). One can see in these expressions for V^ iag the same terms as in (g) 
with exchange — > pj except for the spin-orbit term proportional to j^. 
Before discussing two different strategies for obtaining this last term, let us 
look at the general structure of (§t])- It has the described above matrix form 
and depends on einbein fields p,, i = 1, 2, 3. The latter have been defined 
previously in [ffD[ [IT], as scalars, 2/ij = , and are assumed 

to be found from the stationary point equation in the path-integral form of 
the meson Green's function, or from the stationary point of the Hamiltonian. 
Now the spin-independent part of Hamiltonian is a unit matrix and hence 
can produce scalar stationary values for pj. The situation changes however 
if one tries to incorporate also the SD part of Hamiltonian in the stationary 
point equation for /i, since it would require /1 to have a matrix form similar 
to that of Vsd- 

This is possible in the generalized form of the FFFR, which is now under 
investigation, but in the present form the only possible way of treatment the 
SD part of Hamiltonian is to consider it as a perturbation. For light quarks 
it is not an expansion in l//ij, and the whole expression (^) is obtained 
with the only and numerically good approximation - keeping the bilocal 
(Gaussian) correlator, neglecting all higher ones. 



As it was shown in meson case |20. 28], this perturbation procedure works 



well even for lowest mesons, where SD corrections produce up to around 
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15% of the total mass (a similar situation holds true in earlier quark model 
calculations with fixed and prescribed constituent masses, see e.g. ||). For 
heavier meson and baryon states the masses fa grow rapidly with quantum 



numbers |2(J and validity of perturbative treatment of SD terms becomes 
even better established. In what follows we describe a perturbative procedure 
of treating the nondiagonal terms. 

To this end we must remember (as in point B of derivation in section 4) 
that nondiagonal terms are also present in the preexponential factor (m — D) 
in (|l]). Consider the largest nondiagonal term A^l in (|56|), and take for 
simplicity its asymptotic form 



VW--*E^-A n« = RW/i2«>. (62) 
i=i ^ 

One has (omitting index (i) for simplicity) 



( m _ D) eM -V EE T) 3 ^ ; p ; - ^ j (1 + i-«n + ...T). (63) 

Comparing with the leading term, given by the upper left corner, one 
normalizes (m — D) by extracting the factor [m + \x) and thus obtains 

Bl .(63) = (l-^T + ...), V&> = £ ^g^ . (64, 

In the heavy quark limit, /ij ~ mj and (|H]) coincides with the term 
proportional to in (|8]). For light quarks , when fa 3> rrii however one 
has twice as large coefficient in(|B"TD which coincides with the heavy-quark 
expansion with the light quark expression. Thus our total expression for SD 
potential treated as perturbation is a 2 x 2 matrix 



a 



V^ M quarks) = £ Vjtil + (65) 



r(light quarks) _ , T/( e ) 

r(ij) ■ ■ flcTh /rjm /|COh „„J T/( £ ) 



where is given in (|6T|), ([45]), ( |52|) and Vgjj is given in (|64]), with gen- 
eral form obtained by replacing cr — > J-|. Now one can see from these ex- 
pressions that we have a full correspondence between terms in (|8]) and in 
(|65|), where each term in (^) is obtained from the corresponding one in (H) 



by replacement rrii — > fa, except for the term with J^, where one replaces 
2m, 2 -> /ij(/Xi + rrii). 
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7 Concluding remarks 



We have obtained all perturbative and nonperturbative spin-dependent terms 
in the 3q system in the approximation when lowest (bilocal) field correlator 



is retained in Wilson loop. The analogous procedure for mesons in pq , [19 
yielded SD potentials satisfying Gromes relation |^4H , with correct asymp- 
totics at large distances of Thomas precession type. For the 3q system we 
also get this asymptotics for spin-orbit terms in the form of a sum of one- 
body Thomas terms, in agreement with earlier results in ||26|| . All other 
nonperturbative terms and exact nonasymptotic form of Thomas terms are 
new. The signs of perturbative and nonperturbative spin-orbit terms are 
different and one may expect some cancellation, which should be checked in 
exact calculations of baryon spectra with spin splittings. All nonperturbative 
SD terms in (|8|) except for V5 have the structure similar to that of the QQ 



case, considered in |Tj|, except that spin-orbit terms are of one-body rather 
than the two-body character. The new term V5 (|24| ) does not have a QQ ana- 
log,and after averaging over coordinates has a structure similar to two-body 
spin-orbit force. 

The large N c structure of SD interaction can be clearly seen from explicit 
expressions and may be represented as leading (0(N®)) terms of one-body 
spin-orbit interaction, when both fields in the field correlator are on the same 
sheet of the 3-sheet surface, and suppressed (O^N' 1 )) terms of spin-spin 
interactions and spin-orbit from two different sheets. Hence the 3q dynamics 
in the large N c limit reduces to the uncorrelated motion of N c quarks around 
a common center (string junction), which can be taken as infinitely heavy. 

The general structure of the SD potential (|8|) at large N c is in agreement 
with the classification done in [^UJ where the unsuppressed at large N c terms 
are one-body spin-orbit potentials,while two-body spin-dependent terms are 
1/N C suppressed. In addition in [[^| appear also spin-flavour terms which can 
be associated with with pion and kaon exchange forces. The latter were 
not considered in the present paper, but can be easily included in the same 
formalism, using the new chiral Lagrangian derived in [|30|.It is shown there 
that in the q — string — q system pions are emitted by quarks with the known 
amplitude, so that the pion-exchange force can be predicted unambiguosly 
and added to the those obtained in the present work. This would complete 
the overall picture of SD forces in baryon. The author is grateful to J.Goity 
for useful discussion, remarks and suggestions. 

This work was supported by DOE contract DE-ACOS-84ER 40150 under 



18 



which SURA operates the Thomas Jefferson National Accelerator Facility. 

Appendix A 

Field correlators 



^From general definitions of g 2 (F llv F\ a ) through D,Di in [11] one gets 
q 2 

—tr(F^(x)<$>(x,y)F pa (y)<$>(y,x)) = (5^5 VU - 5 tl J vp )D(z)+ 

+ hd /J ,z p 5 up +perm.]D 1 (z) (A.l) 

a 2 dD 
^-tr(B i (x)^(x,y)B j (y)^(y,x)) = (5 tJ (D(z) + D 1 (z) + z 2 ^]- 

dD 1 

dz 2 

a 2 dD 

j^tr(E i (xMx,y)E j (y)^(y,x)) = (S tJ (D(z) + D 1 (z) + 4~q£] + 

a 2 dD 
jj-tr{B i (x)$(x,y)E j (y)<S>(y,x)) = e ijk z 4 z k — (A.4) 

where we have defined 



(A.2) 



Vn, // = 1,2,3,4, $(x,y) = Pexpig / A^(u)du^ 
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